Expansion of matter waves in static and driven periodic potentials 
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We study the non-equilibrium dynamics of cold atoms held in an optical lattice subjected to a 
periodic driving potential. The expansion of an initially confined atom cloud occurs in two phases: 
an initial quadratic expansion followed by a ballistic behaviour at long times. Accounting for this 
gives a good description of recent experimental results, and provides a robust method to extract the 
effective intersite tunneling from time-of-flight measurements. 
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I. INTRODUCTION 

Experimental progress in confining ultracold atoms in 
optical lattices has undergone spectacular progress in re- 
cent years. Optical lattice potentials are extremely clean 
and controllable, and the excellent coherence properties 
of atomic condensates hold out the prospect of controlling 
their dynamics using quantum coherent methods. This 
is both interesting from the point of view of fundamental 
physics, and has many potential applications to quantum 
information processing, where it is vital that the coher- 
ence of the system is preserved during its time evolution. 
One such scheme is to use a time-periodic driving poten- 
tial to induce the effect termed "dynamical localization" 
[l[ or "coherent destruction of tunneling" This is a 
quantum interference effect in which a particle acquires 
a phase from its interaction with the driving potential, 
which leads to a renormalization of the single-particle 
tunneling probability. For specific values of the driv- 
ing parameters the effective tunneling probability can be 
highly suppressed, providing a sensitive means of coher- 
ently controlling the localization of the atoms [|| . This 
renormalization has recently been directly observed in 
cold atom experiments 041] • 

A convenient way of measuring the effective tunnel- 
ing is to observe the rate of expansion of a conden- 
sate once a harmonic potential trapping the atoms along 
the direction of the optical lattice has been switched off 
d, [j| . Although the results of these experiments agreed 
well with the theoretically expected scaling of the renor- 
malized tunneling probability with the Bessel function 
of the driving strength, for particular initial conditions 
the scaling seemed to be quadratic rather than linear in 
the Bessel function. A number of explanations of this 
phenomenon have since been put forward, including a 
possible crossover from coherent to sequential tunneling 
p, 0] induced by phase scrambling arising from dynam- 
ical instabilities, a time-averaging effect produced by fi- 
nite time- resolution of the measurement [8[ , and driving- 
induced atom-pairing A recent theoretical stability 
analysis (Toj has shown, however, that phase scrambling 



is unlikely to occur for the experimental parameters of 
Refs. [H and the pairing mechanism would require 
rather stronger interactions than were present in those 
experiments. In this paper we propose an alternative ex- 
planation for the observed quadratic scaling of the renor- 
malized tunneling probability based on the exact form of 
the expansion of the condensate, which is linear in the 
long time limit but quadratic for short times fTTj j . Using 
this complete time dependence in order to extract the 
tunneling probability from the experimental data gives 
quantitatively accurate agreement with theory, with no 
adjustable parameters. 



II. MODEL AND ANALYSIS 

The Bose-Hubbard model is described by the Hamil- 
tonian 

#BH = -J J2 [ a i a i + H - c ] + \ J2 n ^ n i - x ) 



where Oj/at are annihilation / creation operators for a 
boson on lattice site j, J (taken to be positive) describes 
the hopping amplitude between nearest neighbor sites 
and U is the repulsive energy between two bosons 
occupying the same site. The operator rij = ajaj is the 
standard number operator, and V(r) is the external trap 
potential, which is usually considered to be parabolic, 
V(r) = mwjr 2 /2, where lot is the trap frequency. Al- 
though simple in appearance, the Bosc-Hubbard model 
can provide an excellent description [lj] of ultracold 
atoms held in optical lattice potentials. 

Adding a static and a sinusoidally varying force to the 
system leads to the general time-dependent potential 



H(t) = H Bn + n i J ( A + K cos ut) . 



(2) 



2 



where A is the static tilt applied to the lattice, and K and 
u) are the amplitude and frequency, respectively, of the os- 
cillating component. Experimentally, the two forces are 
introduced into the rest frame of the optical lattice by ap- 
plying appropriate frequency differences to the acousto- 
optic modulators creating the lattice beams, as described 
in detail in A time-periodic system of this type can 
be analyzed using Floquet theory, revealing 13[ that the 
effect of the driving can be described by the static Hamil- 
tonian (JIJ with a renormalized tunneling J e ff. For an un- 
tilted lattice (A = 0) this renormalization is given by 
the zeroth-order Bessel function J e g = JJq(Kq), where 
for convenience we define K$ = K/hco. Thus at the val- 
ues Ko — 2.404, 5.52, . . ., at which the Bessel function 
vanishes, the effective tunneling is suppressed. This effect 
thus provides a means to coherently control the dynamics 
of trapped atoms, without altering any of the parameters 
of the optical lattice. 

When the trap potential along the lattice direction is 
removed the atom cloud will expand in time, at a rate 
determined by J e ff. To quantify this process we calculate 
the spread of the wavefunction 



a(t) = v^ 2 ) " W 



(3) 



aligning the optical lattice with the x-axis. We begin by 
setting the Hubbard interaction, U, to zero. In the con- 
tinuum approximation, valid when the kinetic energy of 
the condensate is much less than the width of the first 
Bloch band, the ground state of a parabolic trap is simply 
given by a Gaussian, ip(x) = N exp[— x 2 / (2a 2 )], where N 
is the normalization such that /_ \ip(x)\ 2 dx = 1, and 

a = y/W/JpwJr) is the harmonic trap-length. If the trap 
potential is now removed, this initial state will expand 
while remaining Gaussian. This expansion can be calcu- 
lated analytically [11] yielding the result 



<r(t) = a ^Jl + 4{Jt/h) 2 (d L /a) 4 , 



(4) 



where dh is the spacing of the optical lattice, and 
(j o = a/y/2. Similar, but more complicated, expressions 
were obtained by Korsch [141 ] using a lattice represen- 
tation instead of the continuum approximation, which 
coincide with this result. The expansion clearly occurs 
in two different phases, separated by a crossover time 
t c = (h/2\J\) (a/c?i) 2 . For long expansion times, t ^> t c , 
the wavefunction spreads linearly with time, a(t) oc \J\t, 
reproducing the expected ballistic expansion of a released 
wavepacket. For short times, however, t < t c , the expan- 
sion is instead quadratic, a{t) oc J 2 t 2 . It is important 
to note that the expansion depends on both the magni- 
tude of the tunneling, J, as well as on the size of the 
initial wavepacket a. In particular, a tightly-confined 
wavepacket will have a higher spread in momentum, and 
so will enter the regime of linear expansion more quickly. 
The extreme case where only a single lattice site is filled 
was considered in Ref. l], where it was found that the 
expansion was always linear with time. This result is ex- 



actly reproduced by Eq. (j4]) by appropriately taking the 
limit a — > 0. 



III. RESULTS 

Undriven lattice - We first consider the case of a 
static lattice, with A, Kq = 0. In Fig. [T] we show the 
time-dependence of the expansion of an initial Gaussian 
wavepacket, numerically evolved in time under the time- 
dependent Hamiltonian Using J = 0.1E rec , where 
the recoil energy E lec = Ti 2 ir 2 /2md 2 L , we see that the an- 
alytic expression (|4| accords exactly with the numerical 
result, indicating the validity of the continuum approxi- 
mation. The transition from the initial quadratic expan- 
sion to the ballistic regime is clearly visible. Halving the 
value of J used produces the expected result of reducing 
the rate of expansion, and moving the crossover from the 
quadratic to the ballistic regime to a later time. 



Untilted driven lattice - We now consider the ef- 
fect of including the time-dependent driving potential 
V(t) = Kcosut. We choose a high driving frequency, 
fouj = 4J, and tune the amplitude of the driving so 
that </ e ff = 0.5J. In accordance with the predictions 
of the Floquet analysis, we see that the expansion of the 
wavepacket in the driven system with a bare tunneling of 
J = O.lE'rec closely follows the result for J = 0.05-E re c, 
indicating that the driving field indeed renormalizes the 
tunneling as expected. Looking in detail at the expan- 
sion (inset of Fig. [T]) we see that although on average the 
expansion of the driven condensate closely matches that 
of the static case with J = 0. 05-E r ecj the driven result 
contains small amplitude oscillations with the same fre- 
quency as the driving. These oscillations arise from the 
intrinsic time dependence of the Floquet states them- 
selves. Their amplitude reduces as the driving frequency 
becomes larger, indicating that the approximation of 
modelling the driven system with a renormalized static 
Hamiltonian becomes increasingly good. Finally we also 
show in Fig. [1] the most dramatic effect of the renor- 
malization of tunneling. Since J e g = JjTo(-Ko), tun- 
ing Ko to a zero of the Bessel function should result in 
the complete suppression of tunneling (neglecting next- 
nearest neighbor tunneling [Hj]). We indeed see that set- 
ting Ko — 2.404 results in the condensate not expanding 
with time, due to the vanishing of J e ff. Similarly to the 
Ko = 1.22 case, this curve again displays small oscilla- 
tions, which become larger at low values of u>. This low 
frequency behaviour would correspond to the "dynami- 
cal localization" regime [l|, where the wavepacket peri- 
odically returns to its initial state at stroboscopic times 
t = nT = n2n/uj, but between these times can exhibit 
large excursions. 

In Ref. i the effective tunneling was deduced by mea- 
suring the expansion rate of the condensate at a fixed 
time, and assuming this rate was directly proportional to 



FIG. 1. (Color online) Above: Expansion of an initial Gaus- 
sian wavepacket in a flat lattice potential, obtained by the nu- 
merical propagation under Hamiltonian @ . When no driving 
potential is applied (K — 0) a(t) increases following Eq. ([4]). 
For J = O.I-Erec (black solid line) the expansion is clearly 
quadratic initially, and becomes linear at long times. Setting 
J = 0.05i5rec (red solid line) reduces the expansion rate as 
expected. By applying a periodic driving potential the tun- 
neling can be renormalized to an effective value J e ff . Tuning 
Ko = f .22 reduces J e g so that the expansion of the con- 
densate (dashed red line) reproduces the J = 0.05i?rec result. 
Setting Ko = 2.404 - the first zero of the Bessel function - pro- 
duces coherent destruction of tunneling (CDT), and the con- 
densate no longer expands with time (blue dash-dotted line). 
Inset - Detail of the periodically-driven result (Ko = 1.22). 
The driven expansion on average reproduces the J = 0.05-Erec 
result, but shows small oscillations with the same frequency 
of the driving. The amplitude of these oscillations decreases 
with increasing driving frequency. Below: Experimental com- 
parison of the free expansion of a condensate (Ko = 0) with a 
condensate experiencing CDT (Ko = 2.4). As predicted, the 
expansion of the second condensate is strongly suppressed. 



J e ff . Accordingly the ratio between the tunneling param- 
eters in the static and the driven lattice was calculated 
as \ J e s/J\ = (o"(t) - er )/(<7stat - (To), where <7 s tat is the 
size of the condensate after expansion in the static lat- 
tice. For the experimental parameters (d^ = 426 nm and 
J/h = 270 Hz), we calculate a crossover time of t c ~ 9.7 
ms for a weakly-interacting condensate released from a 
20 Hz harmonic trap. As the experiment employed an 
expansion time of 100 ms, the results can thus be ex- 
pected to be reliable only as long as \J e ff/J\ > 0-1- I n 
order to get better agreement with theory we now use 
Eq. containing the full expansion dynamics, giving 

\ JeS /j\ = JW^4- (5) 

V a stat °0 



FIG. 2. Dynamical suppression of tunneling in a periodically 
driven lattice. The effective tunneling parameter was calcu- 
lated assuming a simple linear expansion (open symbols) and 
by taking into account the exact expansion dynamics (filled 
symbols). In the latter case the agreement with the theo- 
retically expected scaling (solid line) is clearly better. The 
experimental parameters were J/h — 240 Hz, uj/2tt = 4kHz 
and the expansion time 150 ms. 

Figure [2] shows that, as expected, using Eq. ([5]) to cal- 
culate the renormalized tunneling gives better agreement 
with theory. 

It is interesting to note that although J e g- is strongly 
suppressed near Kq = 2.4, it does not actually reach zero 
when the Bessel function vanishes. This is due to the 
effect of higher-order hopping terms present in the sys- 
tem's dynamics. Although they too are renormalized by 
the driving potential, for sinusoidal driving they will not 
vanish at the same driving parameters as for the nearest 
neighbor hopping. The residual value of J e g, visible in 
Fig. [21 is in reasonable agreement with the value of the 
next-to-nearest neighbor hopping (around 5%) calculated 
for a similar system in Ref.[15]. 

Tilted lattice, resonant driving - Applying Eq. ^ to 
the experimental data on photon-assisted tunneling Q 
leads to an even more striking improvement in the agree- 
ment between theory and experiment. In those experi- 
ments a tilt was applied to the lattice through a constant 
acceleration, leading to a suppression of tunneling by 
Wannier-Stark localization. Periodic driving of the lat- 
tice at a frequency to matching the energy offset between 
two adjacent lattice wells then led to partial restoration 
of the tunneling probability, with the effective tunneling 
probability given by \ J e g/J\ = Ji{Kq), i-e. one expects 
a scaling with the first-order Bessel function. As shown 
in the inset of Fig. [31 assuming linear expansion in or- 
der to extract J e g/J led to a scaling that interpolated 
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FIG. 3. a) Effective tunneling for resonant driving in a tilted 
lattice (photon-assisted tunneling). The effective tunneling 
parameters were calculated for two different initial condensate 
sizes (around f 5 fim (open symbols) and around 17 /im (filled 
symbols)) using Eq. Q. For comparison, b) shows the same 
experimental data with the renormalized tunneling parameter 
calculated assuming linear expansion. One clearly sees that 
in this case the experimental data interpolate between a lin- 
ear Bessel scaling (solid line) and a quadratic scaling (dashed 
line). 



both data sets give the same dependence on Kq which is 
very close to the theoretical prediction. 



IV. CONCLUSIONS 

We have shown that in order to extract the effec- 
tive tunneling from expansion measurements it is im- 
portant to account for the detailed time-dependence of 
the condensate expansion. When the effective tunneling 
is small, or the initial width of the condensate is large, 
the crossover to ballistic expansion will not be reached 
until very long times. Measurements made at earlier 
times will thus underestimate the effective tunneling rate, 
which gives a quantitatively accurate interpretation of 
the "squared Bessel function" behaviour noted in Refs. 
and@. Although we have not included the effects of in- 
teractions, this is a reasonable approximation for the sys- 
tems studied in Refs. land | where the interactions were 
fairly small (U/h ~ fO Hz), and their effect rapidly be- 
came negligible as the condensate expanded and became 
more dilute. Using the correct expansion formula, Eq. 
(UJ , not only provides an accurate means of deducing the 
value of the effective tunneling, but is essential to study 
subtle effects, such as the transition to diffusive tunneling 
and the influence of higher-order tunneling terms, which 
would otherwise be masked by this behaviour when the 
effective tunneling is small. 
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